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We briefly review how precise determinations of the strong couphng constant and of the heavy 
quark masses may be obtained from heavy quarkonium. Such determinations are competitive with 
heavy quark masses extraction from other systems and give an accurate value for the strong coupling 
constant at a relatively low energy scale. In particular we report about a recent determination of 
cts from r(T(15') — > 7X)/r(T(15') — » X) with CLEO data which includes color octet contributions 
and avoids model dependence in the extraction. The obtained value is as(Af-r(is)) = 0.184lgQj^4 , 
which corresponds to Qs(Mz) = O.IIQIq oos- Future prospects for more precise extractions of the 
heavy quark masses and Qs are discussed. 



I. QCD AND THE ONIA 

QCD is the theory of the strong interactions: we 
should be able to predict all the properties of hadrons 
starting from the QCD Lagrangian which is a func- 
tion only of the coupling constant cks and of the quark 
masses m. Therefore, from the theoretical predictions 
of physical observables and the corresponding exper- 
imental measurements, we should be able to extract 
the values of the coupling constant as and of the quark 
masses m. However, everything is complicated by the 
fact that QCD is a strongly coupled theory in the 
low energy region. At the scale Aqcd nonperturba- 
tive effects become dominant and as becomes large. 
The nonperturbative QCD dynamics originates the 
confinement of quarks that in turn is the reason for 
which the quark mass loses its most intuitive defini- 
tion. Quarks are confined inside hadrons and thus we 
cannot directly measure their masses. The mass of the 
quark is a parameter defined in some renormalisation 
scheme at some renormalisation scale. Systems made 
by two heavy quarks, -quarkonia in the following-, are 
characterized by a quark mass scale mg which is large, 
bigger than Aqcd- Then as(mQ) is small and per- 
turbative expansions may be performed at this scale. 
This introduces a great simplification and hints at 
a factorization between high and low energy contri- 
butions for quarkonia. For these systems, however, 
things are even more interesting [l| . They are nonrela- 
tivistic systems characterized by another small param- 
eter, the heavy-quark velocity u, and by a hierarchy 
of energy scales: mg (hard), the relative momentum 
p ~ rnqv (soft), and the binding energy E ~ rnqv^ 
(ultrasoft). For energy scales close to Aqcd, pertur- 
bation theory breaks down and one has to rely on non- 
perturbative methods. Regardless of this, the nonrel- 
ativistic hierarchy mq 3> rnqv ^ rngv^ will persist 
also below the Aqcd threshold. While the hard scale 
is always larger than Aqcd, different situations may 
arise for the other two scales depending on the con- 
sidered quarkonium system. The soft scale, propor- 
tional to the inverse radius r, may be a perturbative 
(;» Aqcd) or a nonperturbative scale (^ Aqcd) de- 
pending on the physical system. Finally, only for tt 



threshold states the ultrasoft scale may still be per- 
turbative. Heavy quark-antiquark states probe con- 
finement and nonperturbative physics 0] at different 
scales and are thus an ideal and to some extent unique 
laboratory where our understanding of nonperturba- 
tive QCD, its interplay with perturbative QCD and 
the behaviour of the perturbative bound state series 
may be tested and understood in a controlled frame- 
work. In particular in some regimes nonperturbative 
effects will appear in the form of local or nonlocal 
gluon condensates and will be suppressed in the com- 
putation of physical observables. In this framework 
quarkonia become very appropriate systems to be used 
for the study of the transition region from high to low 
energy, for information on the QCD vacuum structure 
and for precision determinations of the QCD param- 
eters. Precisely this last point is the subject of this 
paper. In the next Sections we will discuss the sys- 
tematic framework offered by Non Relativistic Effec- 
tive Field Theories (NR EFT) Q] for the description 
of quarkonia and how one can take advantage of the 
accurate EFT calculations to make precise determi- 
nations of the QCD parameters. For some reviews of 
NR EFTs see 0, i, [ 



II. EFFECTIVE FIELD THEORIES 

It is possible to take advantage from the existence 
of a hierarchy of scales in quarkonia to introduce NR 
EFTs, which are simpler but equivalent to QCD. A 
hierarchy of EFTs may be constructed by systemati- 
cally integrating out modes associated to high energy 
scales not relevant for the quarkonium system. Such 
integration is made in a matching procedure that en- 
forces the complete equivalence between QCD and the 
EFT at a given order of the expansion in v (w^ '--^ 0.1 
for fefe, ~ 0.3 for cc, u - 0.1 for tt). The EFT re- 
alizes a factorization at the Lagrangian level between 
the high energy contributions carried by matching co- 
efficients and the low energy contributions carried by 
the dynamical degrees of freedom. The Poincare sym- 
metry remains intact in a nonlinear realization at the 
level of the NR EFT imposing exact relations among 
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the EFT matching coefficients 0,[i3l- 

By integrating out the hard modes one obtains Non- 
relativistic QCD NRQCD is making exphcit 

at the Lagrangian level the expansions in mv/m and 
mv'^/m. It is is similar to HQET, but with a different 
power counting. It also accounts for contact interac- 
tions between quarks and antiquark pairs (e.g. in de- 
cay processes) and hence has a wider set of operators. 
In NRQCD soft and ultrasoft scales are left dynamical 
and their mixing may complicate calculations, power 
counting and the consideration of the nonperturbative 
effects. In the last few years the problem of system- 
atically treating the remaining dynamical scales in an 
EFT framework has been addressed by several groups 
[ill . IT^ . [l^ and has now reached a good level of un- 
derstanding. Therefore one can go down one step fur- 
ther and integrate out also the soft scale in a match- 
ing procedure to the lowest energy and simplest EFT 
that can be introduced for quarkonia, where only ul- 
trasoft degrees of freedom remain dynamical. We call 
such EFT potential NonRelativistic QCD (pNRQCD) 
[ll[l3i] (an alternative EFT is in [M])- pNRQCD is 
making explicit at the Lagrangian level the expansion 
in mv'^/mv. This EFT is close to a Schrodinger-like 
description of the bound state and hence as simple. 
The bulk of the interaction is carried by potential-like 
terms, but non-potential interactions, associated with 
the propagation of low-energy degrees of freedom {QQ 
colour singlets, QQ colour octets and low energy glu- 
ons), are generally present. They start to contribute 
at NLO in the multipole expansion of the gluon fields 
and are typically related to nonperturbative effects 

In this EFT frame, it is important to establish when 
Aqcd sets in, i.e. when we have to resort to non- 
perturbative methods. For low-lying resonances, it is 
reasonable to assume rnu^ > Aqcd- The system is 
weakly coupled and we may rely on perturbation the- 
ory, for instance, to calculate the potential. In this 
case, we deal with weak coupling pNRQCD. The the- 
oretical challenge here is performing higher-order cal- 
culations and the goal is precision physics. This is the 
case that we will consider in this paper. 

A. The QCD potential and the Static Energy 

The masses may be extracted from a calculation 
of the energy levels and to obtain the energy lev- 
els we need the potential. The QQ potential is a 
Wilson coefficient of pNRQCD [13| obtained by in- 
tegrating out all degrees of freedom but the ultrasoft 
ones. It is given by a series of contributions in an 
expansion in the inverse of the mass of the quark. 
If the quarkonium system is small, the soft scale is 
perturbative and the potentials can be entirely calcu- 
lated in perturbation theory [3|. As matching coeffi- 
cients the potentials undergo renormalization, develop 



a scale dependence and satisfy renormalization group 
equations, which eventually allow to resum potentially 
large logarithms 15]. The static singlet potential (the 
contribution at zero order in the mass expansion) is 
known at three loops apart from the constant term 
[13, [H, HI. The first log related to uhrasoft effects 
arises at three loops. Such logarithm contribution at 
N'^LO and the single logarithm contribution at N^'LO 
may be extracted respectively from a one-loop and 
two-loop calculation in the EFT and have been calcu- 
lated in [13, [HI ■ The static energy given by the sum 
of a constant, the static potential and the ultrasoft 
corrections is free from renormalon ambiguities. By 
comparing it (at the NNLL) with lattice calculations 
one sees that the QCD perturbative series converges 
very nicely to and agrees with the lattice result [26[ 
in the short range and that no nonperturbative linear 
stringy") contribution to the static potential exist 
, [lO]. This is an example of how precise calcula- 
tions may be performed in this framework. Once the 
renormalon contribution has been cancelled, in this 
case between the static potential and the pole mass 
, we are left with a well behaved perturba- 
tive series and we can unambiguosly define power cor- 
rections. It is possible to make predictions of physical 
quantities (in this case the QQ static energy) at high 
order in the perturbative expansion and with a small 
error (including nonperturbative corrections which are 
suppressed in the power counting) and to make a con- 
nection with the lattice results. It is remarkable that 
the dependence on the lattice spacing can be predicted 
in perturbation theory. 



B. The QCD perturbative series of the QQ 
energies and the nonperturbative contributions 

In weak coupling pNRQCD the soft scale is per- 
turbative and the potentials are purely perturbative 
objects. Nonperturbative effects enter energy levels 
and decay calculations in the form of local or nonlo- 
cal electric and magnetic condensates [Hjll^l. We still 
lack a precise and systematic knowledge of such non- 
perturbative purely glue dependent objects. It would 
be important to have for them lattice determinations 
or data extraction (see e.g. [l^) The leading electric 
and magnetic nonlocal correlators may be related to 
the gluelump masses [l^l and to some existing lattice 
(quenched) determinations [1,[13. 

However, since the nonperturbative contributions 
are suppressed in the power counting it is possible to 
obtain good determinations of the masses of the low- 
est quarkonium resonances with purely perturbative 
calculations in the cases in which the perturbative se- 
ries is convergent (after that the appropriate subtrac- 
tions of renormalons have been performed) and large 
logarithms are resummed. In this framework power 
corrections are unambiguously defined. 
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III. nic AND mi, EXTRACTION 

The lowest heavy quarkonium states are suitable 
systems to obtain a precise determination of the mass 
of the heavy quarks b and c. Perturbative determina- 
tions of the T(IS') and J/tp masses have been used to 
extract the b and c masses. These determinations are 
competitive with those coming from different systems 
and different approaches (for the b mass see e.g. [35]). 

Determinations of the quark masses from the per- 
turbative calculation of T and J/V' IS masses differ for 
the order of the perturbative calculation considered, 
for the order of the resummation of the logarithms in v 
and the way in which nonperturbative corrections are 
taken into account. Higher order terms and the resid- 
ual scale dependence of the result give the theoretical 
error on the mass. The main uncertainty in these 
determinations comes from nonperturbative contribu- 
tions (local and nonlocal gluon condensates) together 
with possible effects due to subleading renormalons. 
We report some example of such determinations in 
Tab. 1. 



reference 


order 


fnt{fnb) (GeV) 


[36] 


NNNLO* 


4.210 ±0.090 ±0.025 


[30] 


NNLO +charm 


4.190 ±0.020 ±0.025 


[38] 


NNLO 


4.24 ±0.10 


[37] 


NNNLO* 


4.346 ± 0.070 


[39J 


NNNLO* 


4.20 ± 0.04 


[40] 


NNNLO* 


4.241 ± 0.070 


[41] 


NNLL* 


4.19 ±0.06 


reference 


order 


mcinic) (GeV) 


m 


NNLO 


1.24 ±0.020 


[38] 


NNLO 


1.19 ±0.11 



TABLE L Different recent determinations of mi,(m(,) and 
mc(mc) in the MS scheme from the bottomonium and the 
charmonium systems. The displayed results use either a 
direct calculation of the lowest energy level in perturbation 
theory or non-relativistic sum rules. The * indicates that 
the theoretical input is only partially complete at that 
order. For the detailed discussion about how the error has 
been computed see the original references, for a review see 
i- 

Once the quark masses have been obtained, the 
renormalon subtraction and the same calculational 
approach have been exploited also to obtain the en- 
ergy levels of the lowest resonances. In [2^ a predic- 
tion of the Be mass has been obtained. The NNLO 
calculation with finite charm mass effects [s^ pre- 
dicts a mass of 6307(17) MeV that well matches the 
CDF measurement [3i| and the lattice determination 
[3^ . The same procedure seems to work at NNLO 
even for higher states (inside the theory errors that 
grow) Including logs resummation at NLL, it 



is possible to obtain a prediction for the mass of 
T]t = 9421 ± ll(th)l^((5as) MeV (where the second 
error comes from the uncertainty in ag) and for the 
Be hyperfine separation A = 65 ± 24^^^ MeV [33]. A 
NLO calculation reproduces in part the IP fine split- 
ting m . 

A compilation of values of the b and c mass has 
been presented by the Quarkonium Working Group 
in Chapter 6 of [ij and is reported in Figures [T] and 
[21 The mass determinations presented in such Figures 
include (relativistic and nonrelativistic) sum rule re- 
sults, lattice QCD results, semileptonic B decays as 
well as T(IS') and J/-0 \S determinations. One can 
see that the determinations from quarkonium are com- 
petitive with respect to determinations coming from 
other systems (heavy-light, B decays). The original 
works to which the results in such Figures refer are ex- 
plicitely given and discussed in [l[ . We refer to [l| also 
for an extended review of the different mass schemes, 
the different heavy quark mass extractions approaches 
and the renormalon subtraction. 

From these determinations the QWG reported the 
following values for the MS masses: 

TOb(TOb) = 4.22 ± 0.05 GeV 
mXme) = 1.28 ± 0.05 GeV. 

which are displayed by the darker gray area in Figures 
[1] and O and the following ranges: 

TOfc(rn:h) = 4.12 -4.32 GeV 
mXmc) = 1.18 -1.38 GeV 

corresponding to the lighter gray area in Figures[T]and 
[21 For the details of the calculation of these averages 
and ranges see [l|. 

We see that the QWG values for the b and c mass 
attribute to them an error of 1% and 4% respectively. 
This is a smaller error than the one given in the PDG 
0. 

More recent and more accurate mass determina- 
tions (from lattice unquenched calculation Tnb{rfib) = 
4.4 ± 0.030 GeV [11]; from semileptonic B decays, 
mc(mc) = 1.224 ± 0.017 ± 0.054 GeV [H; from low 
momentum sum rules mf,(mi,) = 4.164 ± 0.025 GeV 
mc(mc) = 1.286 ± 0.013 GeV and a new pre- 
liminary calculation of the mass of the b in the the 
potential subtracted scheme with unquenched lattice 
Fermilab action ^46i] ) would call for a new critical anal- 
ysis and discussion of such extractions and errors and 
an updated mass compilation. 

A. mt from ttbar systems 

In [111, [4^ the total cross section for top quark pair 
production close to threshold in e^e^ annihilation is 
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FIG. 1: Collection of recent bottom quark mass determinations. The circles represent sum rule results, the triangles Up- 
silon IS determinations, the squares lattice QCD results and the upside down triangle a determination from semileptonic 
B decays. The full diamond gives the QWG global average for mi,(mi,). The darker and lighter shaded areas represent 
the QWG error estimates corresponding to a la error and a range respectively. This Table is taken from Chapter 6, pag. 
360 of For a detailed discussion and the explicit references to the original works see 0]. 



investigated at NNLL in the weakly coupled EFT. 
The summation of the large logarithms in the ratio of 
the energy scales significantly reduces the scale depen- 
dence. Studies like these will make feasible a precise 
extractions of the strong coupling, the top mass and 
the top width at a future ILC. The present theoreti- 
cal uncertainties for top mass extraction at the ILC is 
about 100 MeV 



IV. as EXTRACTION FROM QUARKONIA 

The summary of values of as{Mz) from various pro- 
cesses as reported by the PDG 2006 [4^ is given in 
Fig. [31 We see that the value of as as determined from 
quarkonium is considerably smaller than the other de- 
terminations. The effect is seen also in Fig. |4] where 
the values of of as(At) are reported at the values of 
jjL where they are measured. The determination of as 
from T decays is one of the few ones at a relatively low 



energy with a relatively small error. It follows from 
theory calculations of ratio of hadronic and leptonic 
T decays [S^l and use of sum rules for the T system 
[ssl . Isoj , the smaller error being obtained in the first 
case. Here we will report about a determination for 
as from the T decays [64| that has recently solved this 
inconsistency. 

Heavy quarkonium leptonic and non-leptonic inclu- 
sive decay rates have historically provided ways to ex- 
tract as and served as additional confirmation of the 
validity of QCD. Ratios of these quantities are very 
sensitive to as if the data are sufficiently precise. In 
particular, today the inclusive decay widths of J/ip^ 
4'(2S) and T(IS') are known with a few percent error, 
the ones of T(2S'), T(3S') with a 10% error and most 
of the other inclusive decays are known with an error 
of 15-20%. In the last few years the error on char- 
monium P-wave inclusive decays have been reduced 
to half El. On the theory side NRQCD ^ and pN- 
RQCD |49j have provided powerful factorization for- 
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FIG. 2: Collection of recent charm quark mass determinations. The circles represent sum rule results, the triangles 
J/ip IS determinations and the squares quenched lattice QCD results. The full diamond gives the QWG global average 
for mc{fnc)- The darker and lighter shaded areas represent the QWG error estimates corresponding to a la error and 
a range respectively. This Table is taken from Chapter 6, pag. 363 of For a detailed discussion and the explicit 
references to the original works see 



mulas for the inclusive decays. 

S and P wave quarkonium inclusive decays are to- 
day known in the NRQCD factorization up to order 
in the relativistic expansion [1, [HO, HI] and at different 
orders in the perturbative expansion of the matching 
coefficients (see e.g. [i^ for a review). In pNRQCD 
the nonperturbative matrix elements of the four quark 
operators on the quarkonium states can be further 
decomposed in the product of quarkonium wave func- 
tions (or derivatives of quarkonium wave functions) in 
the origin and glue dependent operators, with a sub- 
stantial reduction in the number of nonperturbative 
(and unknown) contributions [i^. A lattice calcula- 
tion of such nonlocal gluonic correlators is however 
still missing. 

Thanks to the EFTs factorization between high en- 
ergy contributions, calculable in QCD perturbation 
theory, and low energy nonperturbative contributions, 
it is possible to consider appropriate ratios of inclu- 
sive decays at some order of the expansion in ag and in 
V. In particular, the ratio T{H ^gg)/T{H ggg) 
(H being a quarkonium state) appears particularly 
promising for the extraction of as ^M, [52, [63] . since 
both the wave function at the origin and the relativis- 
tic corrections cancel out. However, the first mea- 



surements of J/ip and T inclusive radiative decays 
delivered a photon spectrum not compatible with the 
early QCD predictions. Inside the EFT approach it 
was understood that colour octet contributions, ig- 
nored in the early calculations, become very impor- 
tant in the upper end-point region of the spectrum 
[ssj . By considering such octet contributions, using 
pNRQCD to calculate them and Soft CoUinear Ef- 
fective Theory (SCET) to resum end-points singu- 
larities, a good description of the photon spectrum 
has been achieved recently, at least for the T(15') 
state [5J]. These recent theoretical advances com- 
bined with new and more precise data from CLEO 
on T(15) radiative decay j55|, has made the ratio 
= r(T(lS') -> 7X)/r(T(15) ^ X), {X being 
hadrons) particularly suitable for the as extraction at 
the bottom mass scale. For the perturbative calcu- 
lation of the matching coefficients appearing in such 
ratio see [qI, [s^ . Colour octet contributions also affect 
the ratio Rj and are parametrically of the same order 
of the relativistic corrections. They have so far either 
been ignored [s^ or estimated to be small [57]hi the 
available extractions of as from this ratio. In 63], re- 
cent determinations of the T(IS') colour octet matrix 
elements both on the lattice |£0| and in the contin- 
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FIG. 3: Summary of values of Qs(Afz) from various processes, taken from the PDG [43|. The value shown indicate the 
process and the measured value of Qs extrapolated to p = Mz- The error shown is the total error including theoretical 
uncertainties. The PDG average coming from thesee measurements and quoted in the text is also shown. Notice that 
the value of as extracted from T decays is considerably lower than all the other determinations. 



uum [61j have been taken into account. This, together 
with the good theoretical description [54] of the pho- 
ton spectrum (ssj ]. allows for a consistent extraction 
of as{M"[(is)) at NLO. The final result obtained in 
H is 

«s(A/T(i5))=0.184tn}3, (1) 

which corresponds to 

a.(A/z) = 0.119™, (2) 

very close to the central value of the PDG [i^ with 
competitive errors. The key ingredients to get these 
numbers have been the precise CLEO data [5^, the 
use of a QCD calculation to extrapolate the pho- 
ton spectrum at low z, and accurate estimates of the 
colour octet matrix elements, which have been pos- 
sible thanks to recent lattice and continuum calcula- 
tions. At present, the main uncertainty in the extrac- 
tion of as comes from the systematic uncertainties in 

The impact of this determination of as, if included 
in the world average, will be to increase it. 

V. FUTURE PROSPECTS FOR MASS AND 
Qs EXTRACTIONS 

The mass and as determinations from quarkonium 
that we have presented are already competitive with 



the results obtained from other physical systems. 

In the near future as(mc) may be extracted from 
the ratio for the J/tj) provided that a new mea- 
surement of the inclusive photon spectrum for radia- 
tive J/f/' decays will be performed at BESIII [6^. In 
a similar way, the discovery and the measurement of 
the mass with a few MeV accuracy will provide a 
determination of as{Mz) with 3 per mille error from 
the hyperfine separation calculated at NLL [s^ . 

For an improved determination of as from the lat- 
tice calculations of the quarkonium spectrum, we need 
a nonperturbative unquenched determination of Aj^ 
and results on the spectrum obtained with different 
formulations of sea quarks, besides staggered quarks. 
Also the improvement in the lattice extraction of the 
masses would require an improved accuracy in the 
conversion from the bare lattice mass to the MS mass. 
In particular the two loop matching in such conversion 
would be needed for the Fermilab and the NRQCD 
actions. A nonperturbative matching would also be 
desiderable. 

For what concerns the mass extraction from the 
T(IS') and J/ip masses in perturbation theory at 
present, as it has been discussed, the major theo- 
retical error comes from our ignorance of the ultra- 
soft nonperturbative corrections. A lattice calcula- 
tion of the nonperturbative chromoelectric correlator 
together with its matching from lattice to MS scheme 
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FIG. 4: Summary of values of as(/i) at the values of fi where they have been measured, taken from the PDG [43|- The 
line shows the central value and the ilcr limits of the PDG average. The data are in increasing order of jj,: r width, T 
decays, deep inelastic scattering, e^e~ event shapes at 22 GeV from the JADE data, shapes at TRISTAN at 58 GeV, 
Z width, and e'^e~ event shapes at 135 and 189 GeV. Notice how the determination from the T decays is the only one 
outside the band. 



is needed. 

Further improvements in the mass determinations 
from nonrelativistic sum rule would require the full 
NNLL calculation; a complete NNNLO computation 
would also be useful to have a better control on the 
theoretical uncertainties. For low momentum sum 
rules, improved determinations of the R measure- 
ments around bottomonium and charmonium region 
would be crucial. 

We conclude noticing that, within the EFT ap- 
proach and the factorization scheme, precision calcu- 
lations in quarkonium may be applied to all the phys- 
ical observables of the lowest resonances, spectra and 
decays included. To this respect it is particularly in- 
teresting the example of the calculation of Ml transi- 
tions for the lowest quarkonia resonances. In this case 



the Poincare invariance of the EFT imposes exact re- 
lations among matching coefficients that set to zero 
the nonperturbative corrections at order u^. At this 
order the Ml transitions may be exactly calculated in 
perturbation theory [6^ . 
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